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Abstract 

We study the energy spectrum of bottomonium in perturbative QCD, taking a s (Mz) = 
0.1181 ± 0.0020 as input and fixing (m^ ) on the T(1S) mass. Contrary to wide 
beliefs, perturbative QCD reproduces reasonably well the gross structure of the spectrum 
as long as the coupling constant remains smaller than one. We perform a detailed anal- 
ysis and discuss the size of non-perturbative effects. A new qualitative picture on the 
structure of the bottomonium spectrum is provided. The lowest-lying cc and be states 
are also examined. 
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1 Introduction 



In recent years our knowledge of heavy quarkonia in the framework of perturbative QCD has 
developed significantly. On one hand computations of new higher-order corrections to various 
physical quantities appeared |], §, || || . On the other hand the discovery of the mechanism of 
the renormalon cancellation in the quarkonium spectrum || |6| led to a drastic improvement of 
the convergence of the perturbative expansion of the energy levels. As important applications 
up to date, the theory enabled precise determinations of the MS-mass of the bottom quark |7|, |8|, 
H and (in the future) of the top quark [5 1] from (mainly) the energy spectra of the lowest- 
lying states. The main uncertainty comes, in the bottomonium case, from the (essentially) 
unknown non-perturbative contributions. These are generally claimed to be around 100 MeV 
and ultimately set the precision of the prediction. 

As for the structure of the higher quarkonium levels, the Coulomb potential seems to be 
unable to explain it (for a review see []12|). For instance, the spacing between two consecutive 
nS levels does not decrease with n and appears roughly constant. For this reason several 
confining potentials have been introduced in the literature over the last decades. We refer to 



13] for one of the most recent phenomeno logical study on several of them. Progress in this 



direction seems to be no longer possible, as far as the relation of the confining potentials with 
QCD remains obscure. This relation has been elucidated recently in [|14], [L5| where a complete 
and systematic parametrization inside QCD of perturbative and non-perturbative effects of 
the heavy quarkonium interaction has been realized. The parametrization depends on the 
mutual relation between the scale of non-perturbative physics, Aqcd, and the other dynamical 
scales in the specific heavy quarkonium system under study. A way to determine it and the 
size and nature of the non-perturbative contributions, consists in establishing to which extent 
perturbative QCD can consistently and successfully describe the system. This is the aim of 
the present paper, which investigates the range of validity of perturbative QCD on the heavy 
quarkonium spectroscopy and consequently extracts upper bounds to the non-perturbative 
contributions by comparing the perturbative predictions, at the current best accuracy, with 
the experimental data. 

One of the main problems in having a consistent (i.e. convergent) perturbative expansion 
for the quarkonium energy levels has been for a long time its bad behaviour when expressed in 
terms of the pole mass (see, for instance, the poor convergence of the computations reported 
in |16|, 0). The reason has been understood in the renormalon language: the pole mass (T^ 



and the QCD static potential fl8f , respectively, contain renormalon contributions of order 
Aqcd, which get cancelled in the total energy of a color singlet quark-antiquark system || |(|. 
The solution then consists in making explicit this cancellation by substituting in the energy 



expansion a short-range mass (free from infrared ambiguities) for the pole mass. In |L0 



this approach has been used in order to calculate the bottom quark mass by means of sum rules. 
However, an analysis of the consistency between the whole level structure of perturbative QCD 
and the experimental data is lacking and will be done in the present work. In order to make 
explicit the renormalon cancellation we will express the quarkonium energies in terms of the 
MS masses and rearrange the perturbative series in the so-called e expansion Jl9| . This will be 
the key ingredient of our analysis. We note that according to a formal argument, the residual 
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uncertainty of the perturbative expansion due to the next-to-leading renormalon contribution 
is estimated to be Aq CD x (cixAqcd) 2 for a bound-state X of size ax [0 0], e.g. it amounts 
to 5-20 MeV for the IS* bottomonium state. Provided this argument applies, in principle the 
predictions of perturbative QCD can be made precise down to this accuracy by sufficiently 
increasing the orders of the perturbative expansion. 

The paper is organized as follows. In Sec. ^|we set up the formalism and the computational 
strategy. In Sec. |3] we perform the numerical analysis and in Sec. ^ we discuss the errors. In 
Sec. |5] we interpret our result and in Sec. || we draw some conclusions. 



2 Perturbative Expansions 

We assume that the system is non-relativistic, b < 1, where v is the typical size of the 
heavy quark velocity in a quarkonium state. Also, we assume that the soft scale mv is much 
larger than Aqcd- Under these assumptions, the energy (mass) of a quarkonium state X 
is given as a series expansion in the MS coupling constant a s (/x) defined in the theory with 
rii massless quarks only. We write a general expression valid also when the masses of the 
quark and antiquark are different. The reduced mass is defined from the pole masses by 
m r = mi iPole m 2 ,poi e /(mi iPO i e + m 2jPO i e ). Moreover, we define x = 1 - 4m r /(m 1|Pole + m 2jPO i e ); 
when the two masses are equal, x = 0. Up to 0{a^m) the energy of a heavy quarkonium state 
X, identified by the quantum numbers n, I, s and j, is given by* 

E x (fi,a a (fi),m i , 1? d S ) = m ljPO i e + m 2 , P oi e + E U n,x(fi, a s (fj), m iiPO i e ), (1) 

8 2 / ( ) \ ^ 

Eun,x(l^, 0£ s (ji), m^poie) = - ^ a s (/i) 2 m r Y £ k+1 y—^—J Pk(L n i), (2) 

where e — 1 is the parameter that will be used in order to properly organize the perturbative 

expansion in view of the O (Aqcd) renormalon cancellation. Pk(L n i) is a fc-th-degree polynomial 

5 

of L n i = log[3n/i/ (8a s (yu)m r ) ) + S±(n + 1) + -, and the harmonic sums are defined as S p (q) = 
q 1 

2J — • It is convenient to decompose the polynomials into renormalization-group invariant 
subsets: 

= 1, (3) 

Pi = PoLrU + d, (4) 

P2 = ^ 2 ^+(-^o 2 + ^i + ^oCi)L n , + c 2 . (5) 

Here, denote the coefficients of the QCD beta function, (3 = 11 — 2^/3, (3\ = 102 — 38^/3, 
and Cfc's are given by 

* = -4, (6) 



* The full formula up to O(a^m) for the S state spectrum was derived in and later confirmed in Q; 
additional corrections necessary for the spectrum of I > 1 states can be found in Jig] ; the formula for the 
unequal mass case was derived in |M . 



3 



lGvr 2 {2s(s + 1)(1 - x) + 3x} 8vr 2 (D S + 3X L5 ) / . 
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7T 



= g [7(7 + 1) - + 1) - + 1)1 , 

1 1 n 1 n (n + Q! " (n-Z + fc-1)! 

2 21 + J + 2 (n - I - 1)! ^ (n + Z + fc)! A; 3 



(n - Z - 1)! (2Z + fc)! (2fc + 2Z - n) 
+ 2(n + Z)! ^ (A;-l)!(£; + Z-n) 3 ' 



(i 



(7) 

(8) 
(9) 



(10) 



It is understood that the last term of Eq. (^) is zero if n — I < 2.|] 

Next we rewrite the series expansion of E x in terms of the MS masses. This is done by 
expressing the pole masses mj )PO i e in terms of the renormalization-group-invariant MS masses 



rrii — fn iMS {m iMS ): 



m 



i,pole 



4 a s {uii 

mi { 1 + -e 

3 7T 



a Ami) 



7T 



di + e 3 



a Ami 



7T 



The coefficient di is given in ||22|| , while the analytic expression of d 2 can be derived from the 
result of P5[ , the renormalization-group evolutions of a s (/i) and m^(fi), and the matching 
condition Note that the counting in e in Eq. (|2|) and Eq. (|TTj) does not reflect the order 
in a s but the wanted renormalon cancellation. One way to understand this is to consider that 
in the sum of the pole-quark masses and the static QCD potential, X^m^poie + Vqcd(^), the 
renormalon cancellation takes place without reordering of power counting in a s |B|, ||. The 
extra power of a s comes in the energy level expansion when the dynamical variable r~ l is 
replaced by the dynamical scale (nlsj \ r~ l | nlsj ) ~ C pa s m r / Moreover, in order to realize 
the renormalon cancellation at each order of the expansion, it is necessary to expand mj ;PO i e 
and -Ebin.x in the same coupling [25|, |19|, |26|, therefore we express a s (rrTi) in ( pTTj) in terms of 



a s {rrii) = a s (/i) < 1 + 



7T Z 



rrii 



71 



P 2 o 



log 2 (^) + |log 

\m,7 8 



m, 



(12) 



^The infinite sum in Eq. ([10|) can be easily evaluated analytically in terms of £3, etc. for given values of n 
and I, e.g. by using Mathematica. 

* Since the soft scale is, in perturbative QCD, the only scale involved in the calculation of the binding 
energy, -Ebin.x, up to 0(a*m), we may absorb the extra a s by reexpressing -Ebin.x in terms of a s (/x) and the 
soft scales /ij, S oft = ce s ((j,)rfii. Then the order counting in a s (/x) coincides with the order counting in e given 
here. 
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Substituting Eqs. ( |i~2"| ) and ( ITT] ) into Eqs. (g) and (|1|) we get an expression for the energy levels 
of the heavy quarkonium states, which depends on /i, a s (/x) and m^, that we can organize as 
an expansion in e up to order e 3 : 

E x (fi,a s (n),mi) = m 1 + m 2 + E i x\/j,,a s (iJ,) } m i )e + E ( x\/J,a s (iJ,),m i )e 2 

+^ 3) (/i,« s (/i),m i ) £ 3 + ... (13) 

Since the counting in e explicitly realizes the order Aqcd renormalon cancellation, and since 
a s and Wii are short-range quantities, the obtained perturbative expansion (p~3| ) is expected to 
show a better convergence with respect to the original expansion (0). 

The obtained quarkonium mass E x depends on the scale fi, due to our incomplete knowledge 
of the perturbative series. We fix the scale \x by demanding stability against variation of the 
scale: 



= 0. (14) 



When we do this, we expect that the convergence properties of the series become optimal, 
and that the scale becomes close to the inverse of the physical size of the bound-state X. If 
the scale fixed by Eq. ( |I~4] ) evidently does not fulfill these expectations, then the theoretical 
predictions obtained in this way will be considered unreliable. We will show that this typically 
happens when the coupling constant becomes bigger than one. 

Based on the formalism developed in ||14|| , in principle three scenarios are possible under the 
assumption d C 1 and mv 3> Aqcd: 1) The energy or ultrasoft scale mv 2 is much larger than 
Aqcd- In this case the potential is entirely perturbative and non-perturbative corrections are 
parameterized by local condensates of the Voloshin-Leutwyler type pofl . 2) The energy scale 



mv 2 is of order Aqcd- The potential is entirely perturbative and non-perturbative corrections 
are parameterized by non-local condensates. 3) The energy scale, mv 2 , is smaller than Aqcd- In 
this case short-range non-perturbative corrections affect the potential. They are parameterized 
by non-local condensates. A perturbative treatment of the energy levels is consistent only if 
these are small compared to the Coulomb potential. Explicit formulas for the non-perturbative 
contributions in all these three cases can be found in |14|. However, their sizes in an actual 



calculation of the heavy quarkonium spectrum are affected by large uncertainties. The best 
known case is 1): The leading non-perturbative contribution is ~ 1.44n 6 (a s F 2 (0))/(m 3 a s (/x) 4 ). 
It runs out of control for n > 1, so the estimate is suitable only for the heavy quarkonium 
ground states. It is quite sensitive to the value of the gluon condensate and of a s . It is also 
numerically sensitive to the replacement of the strong coupling constant in the MS scheme 
in the denominator with the coupling constant in some other scheme (e.g. the V scheme or 
similar) as used by some authors @, 16]. Non-local condensates, which are suitable for the 



situations 2) and 3), are even less known. Moreover, they get entangled with higher-order 
perturbative corrections (starting at order afmloga s [U) and in particular with the (9(Aq CD ) 



renormalon |14j . Due to these uncertainties a direct evaluation of non-perturbative and higher- 
order contributions will not be included in the present investigation. Neither we will distinguish 
among the three different scenarios outlined above. Our approach consists in looking only at 



5 



the perturbative spectrum up to 0(a^m), under the general assumption that mv ^> Aqcd- 
The internal consistency of the perturbative series (i.e. its convergence) and the comparison 
with the experimental data will provide indirect, but more stringent, contraints on the size of 
non-perturbative and higher-order corrections. 



3 Numerical Analyses 

In this section we examine the series expansions of the quarkonium spectra, Eq. ( 
ically. The input value is a^ 5 \Mz) = 0.1181 ||27|| . We evolve the coupling and match it to 
the couplings of the theory with ni = 4 and 3 successively following and solving the 
renormalizat ion-group equation numerically (4- loop running). In this section we do not take 
into account effects of the non-zero charm mass in the bb and be systems. 

Since we expect the ground states of the bb and cc systems to be the states less affected by 
non-perturbative corrections, we fix and m c through the conditions 

E T(ls) (fi x ,a s (nx),m b ) = E^ 1S) = 9.460 GeV, (15) 
Ej/^{nx,Oi s {nx),mc) = Ej% = 3.097 GeV, (16) 
where the experimental values of the vector ground states have been taken from [pTfl . We 



assume, for the moment, that this identification is not affected by non-perturbative corrections. 
From Eqs. (|15|), ([TB]) and (|14|) we determine \ix (see Tab. |l|), and the b and c MS masses: 

m b = mp(mp) = 4.203 GeV, (17) 
m r = mP(mP) = 1.243 GeV. (18) 



These values are in good agreement with recent estimates based on T H |10| and charmonium 
2T| sum rules respectively. 



Using these masses as input and Eq. (|1J), we can calculate the energy levels of other 
observed quarkonium states. These are given in Tab. |l[ The series expansions for the char- 
monium l 1 ^ state and for the I 3 Si, l 3 Pj, 2 3 Si, 2 3 P and 3 3 Si bottomonium states converge 
well. For these states the differences of the theoretical predictions and the experimental data 
are typically 30-70 MeV. Convergence of the series expansions is poor for the 2Pi, 2P 2 and 
AS bottomonium states as well as for the other charmonium states. We consider that the 
theoretical predictions for these levels are unreliable and marked the corresponding theoretical 
predictions with sharps (Jj) in the case of the bottomonium. In the charmonium case unreli- 
able predictions are not displayed. The differences E^ p — Ex are rather large for the states 
with sharps. Notice that for these states the corresponding a s (fix) becomes larger than one, 
indicating a breakdown of the perturbative series. We also computed the mass of the B^So) 
state. The theoretical prediction is consistent with the experimental value, although the ex- 
perimental error is large. Generally, for states, which we consider reliably calculable in the 
above perturbative approach, the scale dependence decreases as we include more terms of the 
perturbative series. For states, whose predictions we consider unreliable, the series becomes 
much more convergent if we would choose a scale different from (typically larger than) fix ■ 
§ We take the matching scales as Tab and fn c , respectively. 
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State X 


rpexp 
&X 


Ex 


E x p -E x 


e\ 


e\ 


E\ 




f \ 

ac B {Hx) 


J/if) 


3.097 


3.097 





0.362 


0.205 


0.043 


1.07 


0.448 




2.980 


3.056 


-0.076 


0.333 


0.195 


0.042 


1.23 


0.399 


T(l 3 5i) 


9.460 


9.460 





0.837 


0.204 


0.013 


2.49 


0.274 


T(1 3 P ) 


9.860 


9.905 


-0.045 


1.38 


0.115 


0.003 


1.18 


0.409 


T(l 3 Pi) 


9.893 


9.904 


-0.011 


1.40 


0.098 


0.002 


1.15 


0.416 


T(1 3 P 2 ) 


9.913 


9.916 


-0.003 


1.42 


0.086 


0.003 


1.13 


0.422 


T(2 3 Si) 


10.023 


9.966 


+0.057 


1.46 


0.093 


0.009 


1.09 


0.433 


T(2 3 P ) 


10.232 


10.268 


-0.036 


2.37 


-0.66 


0.15 


0.693 


0.691 


T(2 J Pi) 


10.255 


10.316" 


—0.061" 


3.97 


—3.56 


1.50 


0.552 1 ' 


1.20 


T(2 3 P 2 ) 


10.268 


10.457 s 


-0.189 s 


4.55 


-5.03 


2.53 


0.537 s 


1.39 


T(3 3 Si) 


10.355 


10.327 


+0.028 


2.34 


-0.583 


0.163 


0.698 


0.684 


T(4 3 ,S 1 ) 


10.580 


11.760 s 


-1.180 s 


5.45 


-6.47 


4.38 


0.527 s 


1.61 


S c (l^o) 


6.4 ±0.4 


6.324 


0.08 + 0.4 


0.668 


0.187 


0.022 


1.64 


0.329 



Table 1: Comparisons of the theoretical predictions of perturbative QCD and the experimental data. We 
used ni — 4 for 66 systems and ni — 3 for cc and be systems. All dimensionful numbers are in GeV unit. 



The value of /j,x for the l 3 Si state of bb is close to the central value (2.575 GeV) used in 
0, while the value for the l 3 Si state of cc is lower than the central value (1.57 GeV) used by 
the same authors. The value of fix used in Tab. [I] for the l 1 ^ state of the be is close to the 
central value (1.6 GeV) used in [[Zl]]. We also notice a remarkable agreement with their value 
of the B c mass. This supports the idea that the scale of minimal sensitivity is, indeed, close 
to the characteristic scale of the system, which may be identified with the inverse of its size. 

Fig. [I] compares the experimental values of the bottomonium spectrum with theoretical 
predictions after eliminating those states, which are not reliably calculable (2 3 P 1; 2 3 P 2 and 
4 3 Si states). If we take an average of the S'-wave and P-wave levels corresponding to each 
principal quantum number n, the theoretical predictions with a s (Mz) = 0.1181 reproduce the 
experimental values fairly well. On the other hand, the predictions for the S-P splittings and 
the fine splittings are smaller than the experimental values. 



4 Error Estimates 

Besides non-perturbative corrections, there are three different kinds of uncertainties in our 
theoretical predictions for the quarkonium spectra, listed below and in Tab. |[ We examine 
them separately. Also these examinations indicate that the theoretical predictions for the 2 3 Pi, 
2 3 P 2 and 4 3 Si bottomonium states are quite unstable, while the prediction for the 2 3 P state 
is at the boundary. 

1) Uncertainty of ctf\Mz)- We analyzed the quarkonium spectra varying the input pa- 
rameter within the range given in 0: af\M z ) = 0.1181 ± 0.0020. The level spaemgs 
become wider for larger a^(Mz), which is consistent with our naive expectation. We 
find that the central value ajp(Mz) = 0.1181 reproduces the whole level structure better 
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Coulomb Experimental Perturbative 

Values QCD 



Figure 1: The bottomonium level structure as given by a pure Coulombic potential, by experiments and by 
the present analysis: the solid and dashed lines represent the S'-states and P-states, respectively. The input 
parameter of the perturbative QCD calculation is ai 5 ^ (Mz) = 0.1181 (see Tab. |l|). We show only those levels 
that we can compute reliably. The Coulomb levels are calculated with m po i e = 5.105 GeV and a s = 0.5752 
such that they reproduce the 15 and 25* levels. 



than af\M z ) = 0.1201 or 0.1161. 

Charm Mass Effects. We have also made an analysis of the bottomonium spectrum 
including finite charm mass effects. Since they deserve a detailed analysis of their own, 
including new calculations, we will report the full results in a separate paper |3(| . Here we 



only summarize some of the qualitative features of the effects and include them as a part 
of the uncertainties of our present analysis. The corrections to the reliable predictions 
turn out to be positive and to become larger for higher states, ranging up to ~ 200 MeV. 
Much of the effects, however, may be reabsorbed by the uncertainties in a^'(Mz) (as 
given above). 

3) Uncertainties from Higher-Order Corrections. We take the maximum value of the follow- 
ing five estimates as an estimate of uncertainties from unknown higher-order corrections 
for each series expansion: (i) The difference between the theoretical predictions computed 
using a s (fj,) as obtained by solving the renormalization-group equation perturbatively at 



4 loops (Eqs. (3) and (11) of Ref. [28]) and numerically at 4 loops (the data of Tab. 
|l]). (Note that the number of energy levels that can be determined in a reliable way is 
larger with the former definition of the running coupling constant. Also in that case, 
reliable predictions turn out to be close to the experimental values.) (ii) The difference 
between the theoretical predictions computed using the 3-loop and the 4-loop (as in Tab. 
|T|) running coupling constants, fixing o$>{Mz) = 0.1181. (iii) The difference between 
the theoretical estimates obtained by fixing fix on the minimum of \E^\ and the results 
of Tab. m, obtained by fixing fix via the condition ([14]). (iv) The contribution ±|.E^| 
from Tab. [Q. (v) For the IS states we consider the 0{a^m) corrections calculated in the 
large-/?o approximation in [j3|. 



^ We note that the sensitivities of the higher levels to a variation of (x£^{Mz) increase by the charm mass 
effects due to the decoupling of the charm quark. 



8 



For comparison, we list more conservative error estimates. These are the variations of 
rn^c and Ex when we fix the scale as twicefj] of the minimal sensitivity scale: /i = 2 fix, 
where fix is determined from Eq. (|i~4)). 







5af\M z ) 


Estimates of higher-order corrections 






= ±0.0020 


(i) 


(ii) 


(iii) 


(iv) 


(v) 


±max 


\i = 2/i X 




-19 
+18 
-16 
+15 


-2 
-5 


+1 
+4 




+3 


±7 
±21 


+5 
+20 


±7 
±21 


+16 
+37 


5E X 


T(1 3 P ) 

T(1 3 P!) 

T(1 3 P 2 ) 
T(2 3 5i) 
T(2 3 P ) 
T(2 3 Pi) 
T(2 3 P 2 ) 
T(3 3 5i) 
T(4 3 Sx) 
5 C (1%) 


+54 
-48 
+63 
-42 
+57 
-50 
+65 
-58 
+117 
+755» 
+83" 
+1189« 
-44" 
+1528" 
+130 
-93 
-381" 
+1883« 
+4 
-5 


+ 15 
+22 
+16 
+18 
+33 
-4 

-136 
+37 

-639* 
-4 


-14 
-8 
-16 
-19 
-57 
+ 1637» 
+2136" 
-63 
+2936" 




+7 

-1 
-21 
-66* 
-206* 
-36 
-1425" 
-1 


±3 
±2 
±3 
±9 
±150 
±1500* 
±2530* 
±163 
±4380* 
±22 




±15 
±22 
±16 
±19 
±150 
±1637» 
±2530* 
±163 
±4380* 
±22 


-53 
-48 
-54 
-73 
-120 
-97 
-229 
-161 
-1361 
+ 1 



Table 2: Variations of the theoretical predictions of Tab. |T] when the uncertainties 1) and 3) discussed in Sec. 
H are separately taken into account. All dimensionful numbers are in MeV unit. Those values corresponding to 
the unreliable theoretical predictions are marked with sharps. The input parameter is taken as a s (Mz) = 0.1181 
except in the first column. The column "±max" lists ±max{|(i)|, |(ii)|, |(iii)|, |(iv)|, |(v)|}. In the last column 
we write conservative estimates with the scale choice fi — 2/ix ■ 



5 Interpretation 

The most non-trivial feature of the present theoretical predictions for the bottomonium spec- 
trum is that the level spacings between consecutive n's are almost constant, whereas in the 
Coulomb spectrum the level spacings decrease as l/n 2 |] Conventionally, this same difference 
between the Coulomb spectrum and the observed quarkonium spectra motivated people to 
construct various potential models. It is, therefore, imperative to elucidate how the above 
perturbative QCD calculation is able to reproduce such a level structure. We will focus on two 
points: (1) the leading renormalon cancellation, which implies that infrared physics decouples; 
this is essential to obtain convergent series expansions; (2) the meaning of the scale fix chosen 
by the scale-fixing prescription flT4]). 

II If we fix the scale as half of the minimal sensitivity scale, fi = (Ax/2, the predictions for the n = 2 
bottomonium states appear to be meaningless, since the scales are quite close to the infrared singularity of the 
runnning coupling constant, and the predictions for the n > 3 states do not exist, since the scales lie below the 
infrared singularity. 

* If we consider, for instance, the ratio {E 3 s — -E^sOA-E^s — Eis), then we obtain experimentally 0.59, from 
the data of Tab. |], 0.71, while from a pure Coulomb interaction 0.19. 
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Figure 2: The support functions fx(q) vs. q ior X — 15, 25 and 35 (solid lines). fx(q) is calculated using 
TOpoic = 5 GeV and a different a s (/ix), taken from Tab. [l], for each X . Vertical lines represent the corresponding 
scales \ix taken from the same table. Also ai 4 \q) is shown by a dashed line. 



Let us approximate each term on the right-hand-side of Eq. (|1|) (in the equal-mass case) as 

d 3 q 



2m 



pole 



2m + 



2CV r w 
|V QC d(?)| = 2m+ ~ 



7T JO 



dqa v (q), 



bin,X 



(X 



\q\<m 
2 



+ V QC d\X). 



(19) 
(20) 



Here, Vqcd(?) = —Cp^Oiv{q) / q 2 is the QCD static potential in momentum space; \X) de- 
notes the Coulomb wave function (the zeroth-order approximation) and Cp = 4/3. The first 
approximation follows from the fact that the dominant contribution to the pole-MS mass rela- 
tion can be read from the infrared region, loop momenta g < m, of the QCD static potential 
0. The second approximation is more obvious. Let us note that also the right-hand-side of 
Eq. (|20D incorporates and is dominated by the leading renormalon contribution to the static 
potential. From Eqs. (|T9|) and ( p0|) we obtain 



Ex ~2m+—^ dqa v {q)f x {q) + (X\^—\X)~2m+^ dq a v (q) f x (q). (21) 
n Jo rn po i c n Jo 

The last approximation follows from the fact that the kinetic energy contribution to the bot- 
tomonium levels turns out to be numerically small[| (notice that this does not contradict the 
virial theorem, since the static potential we are considering here incorporates the running of 
a s and therefore is not simply of the form 1/r). The support function fx is 



fx{q)=0{m-q)- drr 2 \R x (r) 
Jo 



sin(gr) 



qr 



(22) 



where Rx{f) is the radial part of the Coulomb wave-function of X. fx{<l) is almost unity 
in the region l/ax ~ q < fn, where ax is the inverse of the dumping scale of fx and may 

tit is about 17% of for the X = IS state, 6% of for the X = 25 state, 4% of E^ ] for the X = 35 



state. Moreover, these contributions tend to cancel each other in the level spacings. 
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be interpreted as the size of the bound-state X!^\ For the IS" state fx(q) dumps slowly as 
q decreases. For other states fx(l) dumps rapidly from scales which are somewhat smaller 
than the naive expectations {CFa s m po i c )/nx- Eq. ([ZT|) tells that the major contribution to the 
energy levels comes from the region 1 /ax ^ q < m of the self-energy corrections of quark 
and antiquark (apart from the constant contribution 2m). In Fig. |2| we show fx for different 
states calculated with m po i e = 5 GeV and a s (/ix) taken from Tab. |l| The corresponding values 
of nx are also displayed. For those states which we consider the predictions reliable, /ix is 
located within the range where fx(o) — 1 (close to the infrared edge); for those states with 
unreliable predictions, /ix is located out of this range but far in the infrared region. In Fig. Q 
also a s (q) is shown. We see that as nx increases from 1 to 3, the coupling a s (q), close to the 
dumping scale of fxio)-, grows rapidly. According to Eq. (^T[), this is the very reason for the 
widening of the level spacings in excited states in comparison to the Coulomb spectrum. 

Summarizing the above discussion we may draw the following qualitative picture of the 
structure of the bottomonium spectrum: 

• The energy levels of bottomonium are mainly determined by (i) the MS masses of b and 
b, and (ii) contributions to the self-energies of b and b from gluons with wavelengths 
1/m < A < ax- The latter contributions may be regarded as the difference between the 
(state-dependent) constituent quark masses and the current quark masses. 

• Level spacing between consecutive ra's increases rapidly with n as compared to the 
Coulomb spectra. This is because the self-energy contributions (ii) grow rapidly as the 
physical size of the bound-state increases. 

6 Conclusions 

For all the bottomonium states, where the predictions of perturbative QCD can be made 
reliably (i.e. a 8 < 1), our results are consistent with the experimental data within the estimated 
uncertainties of the theoretical predictions. The obtained value for % is in good agreement with 
the recent sum-rule calculations. The theoretical uncertainties given in Tab. [| are numerically 
of the same size as Aqcd X (^xAqcd) 2 , i-e. of the effect of the C(Aq CD ) renormalons: if we 
approximate l/ax — Hx, take the values of Tab. [I], and A QCD = 300 — 500 MeV, we obtain for 
the IS" state a contribution of order ±(5 — 20) MeV, for the n = 2 states a contribution of order 
±(20 - 110) MeV and for the 3S state a contribution of order ±(50 - 250) MeV. Since the 
mass m b has been fixed on the vector ground state and has not been adjusted for higher states, 
the data at our disposal suggest that: 1) the bulk of the bottomonium spectrum is accessible 
by perturbative QCD up to some of the n = 3 states; 2) non-perturbative contributions do not 
need to be larger than 250 MeV for the reliable n = 3 states, than 100 MeV for the n = 2 states 
and than 20 MeV for mi, and may be of the type associated with the (9(Aq CD ) renormalon. 
These upper bounds to the non-perturbative corrections are conservative and their true sizes 

• According to Eq. ( pl| ) , l/ax acts as an infrared cut off in the computation of the energy level. We may 
compare it with a qualitative picture where the gluons, whose wavelengths are much larger than the size of the 
color-singlet bound-state, cannot couple to it. 
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may be considerably smaller; note that for reliable predictions all of \E e ^ v — Ex\ in Tab. |1] 
are smaller than 80 MeV. The existence of non-perturbative corrections of the above sizes may 
also explain the discrepancy observed in Tab. |I| between the theoretical estimates and the 
experimental data of the fine structure and the 25'- IP splittings. However, before making any 
definitive statement, an analysis including higher-order corrections is necessary. 

The reliability of the perturbative calculations and the agreement with the experimental 
data indicates that the assumptions made at the end of Sec. § are, indeed, satisfied by the n — 1, 
n = 2 and some of the n = 3 bottomonium states. Hence, non-perturbative contributions are 
encoded into non-local condensates, which may reduce to local ones for the ground states. These 
are expected to reabsorb the C(Aq CD ) renormalon [TJj] and are of the order Aqcd x (g-xAqcd) 2 - 



The agreement noticed above, between this estimate and the uncertainties of the perturbative 
series given in Tab. 0, shows the consistency of our conclusions. 

For what concerns the cc system, with the present method we are not in the condition to 
make reliable predictions for states higher than the ground state, and, therefore, we cannot 
extrapolate from consistency arguments the size (and the nature) of the non-perturbative 
corrections. We have noticed, however, that our estimate for m c is in good agreement with 
recent sum-rule calculations. From the prediction of the rj c mass we may guess non-perturbative 
contributions to be, in this case, of the order of 100 MeV. Also this figure is consistent with 
the C(A^ CD ) renormalon effect (20 - 110 MeV). 
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